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SUMMARY 


The  theory  for  single-degree  torsional  instability  of  a  two- 
dimensional  wing  oscillating  in  a  supersonic  stream  is  extended  so  as 
to  apply  to  a  finite  rectangular  wing  oscillating  in  a  supersonic  stream. 
The  velocity  potential  and  aerodynamic— torsional-moment  coefficient 
based  on  the  linearized  equations  of  motion  for  small  disturbances  are 
derived  by  means  of  appropriate  distributions  of  moving  sources  and 
doublets.  The  aerodynamic— torsional-moment  coefficient  thus  derived  is 
combined  with  a  mechanical— damping  coefficient  to  study  the  effect  of 
aspect  ratio  on  the  undamped  torsional  oscillations  of  a  finite  rec¬ 
tangular  wing.  Decreasing  the  aspect  ratio  of  the  wing  is  found  to 
have  a  highly  stabilizing  effect  on  the  undamped  torsional  oscillations. 
Results  of  some  selected  calculations  are  presented  in  several  figures. 

It  is  pointed  out  that  second— order  thickness  effects  may  be  of 
significance . 


INTRODUCTION 


In  theoretical  studies  of  an  oscillating  wing  in  two-dimensional 
supersonic  flow  Possio  noted,  reference  1,  that  under  certain  conditions 
a  single— degree  torsional  instability  is  possible.  This  instability, 
also  sometimes  known  as  a  type  of  "single-degree  flutter,"  was  briefly 
discussed  by  Temple  and  Jahn  in  reference  2  and  has  since  been  further 
investigated  by  Garrick  and  Rubinow  in  reference  3j  hy  W.  P.  Jones  in 
reference  4,  and  by  Cheilek  and  Frissel  in  reference  5. 

It  is  pointed  out  in  reference  3  that  this  single— degree  flutter 
is  due  to  the  wing  being  negatively  damped  in  torsion  and  that  the 
negative  damping  is  associated  with  a  change  in  sign  of  the  torsional- 
damping  coefficient;  furthermore  in  the  two-dimensional  case  the 
instability  may  take  place  in  a  low  supersonic  Mach  number 
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range  (l  <Mach  number  |/275),  at  low  values  of  the  frequency  parameter, 
and  for  axis-of— rotation  locations  forward  of  a  point  two— thirds  of  the 
chord  distance  from  the  leading  edge.  This  phenomenon  may  he  of 
particular  importance  in  connection  with  high-speed  airplanes  that  are 
flown  in  (or  through)  the  range  of  low  supersonic  Mach  numbers  at  which 
the  phenomenon  may  take  place.  Although  in  general  discussions  this 
phenomenon  is  usually  associated  with  the  word  "wing*1  the  theory  applies 
as  well  to  an  aileron  when  the  aileron  is  considered  as  a  separate 
degree  of  freedom.  Hereafter  in  this  discussion  the  word  "wing"  may 
generally  he  given  the  broader  interpretation  of  "wing  or  aileron." 

The  purpose  of  the  present  paper  is  to  extend  the  theoretical 
investigation  of  the  single-degree  torsional  instability  of  a  wing 
oscillating  in  two-dimensional  supersonic  flow,  as  presented  in 
reference  3,  to  a  finite  rectangular  wing  and  to  determine  the  effect 
of  aspect  ratio  on  the  undamped  torsional  oscillations  of  such  a  wing. 

It  is  assumed  that  the  negative  damping  phenomenon  in  three-dimensional 
flow,  like  that  in  two-dimensional  flow,  is  determined  mainly  by  low- 
order  terms  of  a  low  frequency  and  that  only  the  effect  of  the  first 
power  of  the  frequency  must  be  considered  to  get  a  good  approximation 
of  its  total  effect. 

In  order  to  obtain  the  three-dimensional  velocity  potential  and 
aerodynamic— damping  coefficient  the  method  suggested  by  Garrick  and 
Rubinow  in  reference  6,  which  is  briefly  discussed  in  subsequent 
paragraphs,  is  applied  to  a  thin,  flat,  rectangular  wing  performing 
slow,  sinusoidal,  torsional  oscillations  in  a  supersonic  stream.  The 
particular  wing  treated  is  such  that  the  Mach  cones  emanating  from  the 
foremost  point  of  each  tip  do  not  intersect  the  opposite  tip  ahead  of 
the  trailing  edge  of  the  wing. 

The  procedure  developed  herein  may  be  readily  extended  to  apply  to 
any  plan  form  with  supersonic  trailing  edge  as  long  as  other  edges  that 
might  be  in  the  regions  of  mixed  supersonic  flow  are  continuously  straight. 
What  might  be  a  more  desirable  extension  would  be  to  obtain  the  non¬ 
linear  effect  that  thickness  might  have  on  the  undamped  torsional  oscil¬ 
lations  of  a  given  plan  form. 

In  reference  6,  Garrick  and  Rubinow  make  use  of  the  theory  of  small 
perturbations  to  investigate  the  air  forces  on  a  thin  finite  wing 
oscillating  in  a  supersonic  stream.  For  convenience  the  boundary— value 
problems  for  the  velocity  potential  for  a  three-dimensional  surface 
(finite  wing)  moving  at  supersonic  speed  are  classified  into  two  types 
and  referred  to  as  "purely  supersonic"  and  "mixed  supersonic."  The 
purely  supersonic  boundary— value  problem  refers  to  regions  of  flow 
where  no  interaction  between  the  flow  on  the  upper  and  lower  surfaces  is 
present.  In  this  case  the  surfaces  can  be  treated  separately  and  the 
boundary— value  problem  for  each  surface  can  be  satisfied  by  source  and 
sink  distributions.  The  source  and  sink  distributions  for  each  surface 
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are  known  f unctions  of  the  plan  form  and  profile  of  the  vingo  Conse¬ 
quently,  the  velocity  potentials  in  the  purely  supersonic  region  can 
always  he  expressed  in  the  form  of  surface  integrals  with  known 
integrands „ 

The  mixed-supersonic  boundary— value  problem  refers  to  regions  where 
interaction  between  the  flow  on  the  upper  and  lower  surfaces  is  present,, 
This  problem  cannot  be  solved  by  a  distribution  of  souroes  alone  over 
the  wing  surface  but  in  general  can  be  satisfied  by  use  of  doublet 
distributions o  The  manner  in  which  the  doublets  are  to  be  distributed 
(or  the  distribution  function)  depends  on  the  camber  of  the  wing  and  on 
the  plan  form  of  the  region  of  the  wing  where  interaction  between  the 
flow  on  the  upper  and  lower  surfaces  is  present.  In  order  to  find  a 
required  distribution  function  it  is  usually  possible  to  make  use  of  the 
given  boundary  conditions  and  express  the  distribution  function  as  the 
unknown  function  in  an  integral  equation. 

It  may  be  appropriate  to  mention  that  Evvard,  reference  7,  has 
recently  developed,  by  consideration  of  a  source  distribution  over  the 
entire  upwash  field.,  a  time— dependent  velocity  potential  that  may  be 
applied  to  certain  edge  problems.  In  reference  8  Harmon  made  use  of 
this  development  to  derive  some  stability  derivatives  for  thin 
rectangular  wings  at  supersonic  speeds.  The  velocity  potential  derived 
herein  for  the  slowly  oscillating  case  can  be  shown  to  bear  a  relation¬ 
ship  to  the  sum  of  three  potentials  employed  by  Harmon,  namely  those 
due  to  (l)  constant  vertical  motion,  (2)  accelerated  vertical  motion, 
and  (3)  pitching  motion 0 
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disturbance— velocity  potential 

potential  function  of  a  moving  source  defined  in 
equation  (7) 

potential  function  of  a  moving  doublet  defined  in 
equation  (l6) 

potential  function  due  to  distribution  of  sources  in 
region  N  of  figure  1 

potential  function  due  to  distribution  of  doublets  in 
region  T  of  figure  1 


x,  y,  z  rectangular  coordinates  attached  to  wing  moving  in  nega¬ 

tive  x-direction 
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function  defining  mean  ordinates  of  any  chordwise  section 
of  wing  such  as  at  y  =  y^  as  shown  in  figure  1 


w(x,yx,t) 


vertical  velocity  at  surface 


of  wing  along  chordwise 


section  at 


n 


X 


0 


t 

a 

• 

a 

V 

c 

M 


rectangular  coordinates  used  to  represent  space  location 
of  sources  or  doublets  in  the  xy— plane 

abscissa  of  axis  of  rotation  of  wing  (elastic  axis)  as 
shown  in  figure  1 

time 

angle  of  attack 

time  derivative  of  a 

velocity  of  main  stream 

velocity  of  sound 

free— stream  Mach  number  (V/c) 


p  =  \Ju2  -  1 

Tf,  t ^  flmc',:;ions  defined  with  equation  (7) 


function  used  to  represent  space  variation  of  source  and 
doublet  strengths 


w(t  -  t-,_)  1 
w(t  -  t2)  J 


functions  used  to  represent  time  variation  of  source  and 
doublet  strengths 


CD 


angular  frequency 


b 


one— half  cord 


h 


one— half  span 


4 
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A 


aspect  ratio 


h 
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F(|,t)),  f(|,T])  functions  used  to  denote  doublet  distributions 
p  density 

p  pressure  difference,,  measured  positive  downwards,  defined 

in  equation  (24) 

Pw/p  values  of  p  in  regions  N  and  T,  respectively,  of 

T  figure  1 


m4 

V 


P 

m 

r, 

( 

I 


aerodynamic  moment.,  defined  in  equation  (25) 

in—phase  component  of  aerodynamic-moment  coefficient 
defined  in  equation  (30) 

out— of— phase  component  or  (torsional-damping)  aerodynamic- 
moment  coefficient ^  defined  in  equation  (31) 

total— torsional— damping  coefficient 

mechanical— damping  coefficient;  rtg  corresponds  to  the 
usual  logarithmic  decrement 

wing  density  parameter  ( 

\4b2 

mass  of  wing 

radius  of  gyration  divided  hy  h; 

moment  of  inertia  per  unit  length  of  wing  about  elastic 
axis 


a  =  M 

Y 

Y 

-  flutter  coefficient 

03a  natural  angular  frequency  of  torsional  vibration  about 

elastic  axis 
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ANALYSIS 


Consider  a  thin  rectangular  wing  moving  at  a  constant  supersonic 
speed  in  a  chordwise  direction  normal  to  its  leading  edge.  In  accord¬ 
ance  with  linear  theory  the  boundary— value  problem  for  the  velocity 
potential  is  treated  at  the  plane  of  the  wing.  For  the  portion  of  the 
wing  between  the  Mach  cones  emanating  from  the  foremost  point  of  each 
tip,  region  N  in  figure  1(a),  the  boundary— value  problem^ involves  only 
the  purely  supersonic  flow;  and  for  the  portions  of  the  wing  within  these 
Mach  cones,  regions  T  in  figure  1(a),  the  boundary— value  problem  involves 
the  mixed  supersonic  flow.  Thus  the  solution  to  the  present  problem  pro¬ 
poses  the  use  of  the  two  types  of  boundary— value  problems  discussed  in 
the  introduction. 


Boundary— value  problems  for  velocity  potential  for  a  rectangular 
wing.—  The  differential  equation  for  the  propagation  of  small  disturbances 
that  must  be  satisfied  by  the  velocity  potential  in  both  regions  T 
and  N  (referred  to  a  uniformly  moving  coordinate  system  as  shown  in 
fig.  1)  is  (equation  4,  reference  6) 


1_ 

2 


r-  +  V 


=  £2  +  + 
dx2  dy2  dz2 


(1) 


The  boundary  conditions  that  must  be  satisfied  by  the  velocity  potentials 
may  be  stated  as  follows:  (a)  In  regions  T  and  N  the  flow  must  at 
all  times  be  tangent  to  the  wing  surface  and  (b)  in  regions  T  the  pressure 
must  fall  to  zero  along  the  wing  tips  and  remain  zero  in  the  portion  of 
the  Mach  cones  emanating  from  the  foremost  points  of  the  wing  tips  not 
occupied  by  the  wing.  (As  customary  with  linear  boundary  conditions  the 
effects  of  thickness  and  camber  are  separated.  In  particular,  in  the 
nonstationary  case  the  linear  thickness  effect  is  of  no  significance  and 
is  not  considered  here0  The  important  camber  effect  is  characterized  by 
the  conditions  that  the  perturbation  pressure,  the  perturbation  velocity 
in  the  free— stream  direction,  and  the  perturbation-velocity  potential  are 
all  anti symmetrical  with  respect  to  the  reference  plane,  whereas  the 
perturbation  velocity  normal  to  the  reference  plane  is  symmetrical 0 
Accordingly,  the  boundary  condition  of  zero  pressure  in  the  side  wake  may 
be  stated  as  zero  perturbation  velocity  in  the  free— stream  direction.  In 
view  of  these  antisymmetric  characteristics  It  is  only  necessary  to  derive 
the  perturbation— velocity  potentials  for  either  the  upper  or  lower  surface 
of  the  wing.  The  upper  surface  is  chosen  for  the  derivations  herein.) 

With  the  boundary  conditions  stated  as  in  the  preceding  paragraph, 
the  difference  in  the  two  types  of  boundary— value  problems  involved 
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is  in  the  additional  conditions  stipulated  in  condition  (b)  that  must 
be  satisfied  in  region  To 

The  condition  of  tangential  flow  may  be  expressed  analytically  as 


(|f)  -  wfcTi.t) 


8x  +  8t 


(2) 


where  Zm  is  the  function  defining  the  mean  ordinates  at  any  chordwise 

section  of  the  wingo  For  the  particular  case  of  the  wing  performing 
small  sinusoidal  oscillations  of  maximum  amplitude  aQ  about  the  span- 

wise  axis  x  =  Xqj  the  equation  for  may  be  written  as  (see  fig.  1(b)) 


(3) 


Substituting  this  expression  for  Zm  into  equation  (2)  gives  for  this 
case 


w(xj,t) 


¥a  + 


d(x  -  xQ) 


(4) 


Equation  (4)  indicates  that  the  vertical  motion  of  a  wing  moving  for¬ 
ward  and  at  the  same  time  performing  sinusoidal  oscillations  about  a 
spanwise  axis  is  equivalent  to  the  superposition  of  the  vertical  motion 
of  two  similar  wings  moving  forward.,  one  at  an  instantaneous  angle  of 
attach  and  the  other  rotating  about  the  spanwise  axis  x  =  Xq  at  an 

instantaneous  rate  of  rotation Q  Since  the  differential  equation  (l) 

is  linear,,  the  velocity  potential  satisfying  equations  (l)  and  (4)  may 

be  considered  as  the  sum  of  two  potentials:  the  first  corresponding 

to  the  first  term  of  the  right-hand  side  of  equation  (4)  which 

hereinafter  is  denoted  by  j)  and  the  second  corresponding  to  the 

co 

second  term  of  the  right-hand  side  of  equation  (4)  which  hereinafter 
is  denoted  by  0^.  The  symbolic  forms  of  the  velocity  potentials  for 
the  regions  N  and  T  may  therefore  be  written  as 


0=0 


Na 


(5) 
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and 

=  ^Ta  +  ^Ta 


A 


(6) 


Derivation  of  0^.—  Reference  6  shows  that  the  "boundary— value 

problems  for  unsteady  motion  in  the  purely  supersonic  region  (region  W) 
may  he  satisfied  by  distributions  of  moving  sources.  The  potential  at 
any  point  (x^y,  z)  due  to  a  moving  source  of  strength  w(t)J,  varying  only 
with  time,,  located  at  point  (^r],0)  is  given  in  equation  (7a)  of  refer¬ 
ence  6.  In  the  present  notation  the  expression  for  this  potential  may 
be  written  as 

K  -  (T) 

\Kn  -  -  ") 

where 

t  =  M(x  -  j )  _  ^ )  (n?  -  n) 

1  cp2  pc  « 


t g  =  Mix  -  §.)  + 
cp2 


\  =  y  -  |  \/(x  -  l  )2  -  p2z2 

112  =  y  +  i  V&  -  5  )2  -  P2z2 


* 
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The  integral  form  of  the  velocity  potential  at  any  point  (x,y,z )  due  to 
a  distribution  of  such  sources  over  the  |  tj- plane  is 


ix-pz  nr]2 

W(|,tj)0i  dll  d£  (8) 

3  Jni 

where  W(|,i})  represents  the  space  variation  of  source  strength. 

The  limits  tj  and  r\^  are  the  intersections  of  the  |  tj— plane 

and  the  upstream-opening  Mach  cone  with  vertex  at  point  (x,y, z).  Thus 
the  velocity  potential  at  point  (x,y,z)  is  affected  only  hy  the  sources 
in  the  |T)-plane  that  are  within  this  Mach  cone.  (See  fig.  2(a)  for 
the  limiting  case  z  =0.) 

In  order  to  obtain  for  an  oscillating  flat  wing,  the  time 

variation  of  source  strength  w(t)  is  defined,  as  follows: 

w(t)  =  e1(Dt 


0(x,y,z,t)  =  - 


2  pit 


(i 


This  form  for  w(t)  gives  for  the  numerator  of  equation  (7) 


w(t  -  T3_)  +  w(t  —  T 2) 


ia)(t— ' Ti)  iao(t— tp) 

e  N  +  e  N  ~ 


T2+Tl 

•  10j  To 

~  2ela^  e  2  cos  a>- ■ — — -  (9) 

2 


The  space  variation  of  source  strength  is  defined  as 


W  (l.ri)  =  -Vcc0 


(10) 
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If  the  expressions  given  in  equations  (9)  and  (10)  are  substituted  into 
equation  (8)  and  z  is.  made  zero,  the  following  equation  for  the  velocity 

potential  on  the  upper  surface  of  the  wing,  with  a  =  a0eiayfc,  is  obtained: 


ioM 


0Na(x^^)  =  7? 


(x-|) 


COS 


W  -  ii)(ig  -  n) 


\/(i  -  la)  (\  - i) 


dr)  d£ 


(ID 


4 


The  integrations  indicated  in  equation  (ll)  can  be  carried  out  in 
the  form  of  a  series  of  Bessel  functions  as  in  reference  3;  however,  this 
is  not  required  here  since,  as  is  pointed  out  in  the  introduction,  only  the 
first-order  effects  of  the  frequency  cd  are  necessary  iii  the  present 
discussion.  The  terms  involving  higher  powers  of  the  frequency  can  be 
deleted  by  expanding  the  integrand  of  equation  ^11)  into  a  power  series 
m  co  and  dropping  the  terms  involving  this  parameter  to  powers  higher 
than  the  first.  When  this  expansion  is  made,  the  indicated  integration 
m  equation  (11)  can  be  carried  out  in  closed  form.  There  is  obtained 


0Na 


<*-t) 

CB2 

^  -  \)("2  -  ") 


i  dr]  d  I 


Vcc 

P 


IcM 

2P2c 


I 


ori  "the  out/ — of* — pha s e  component 


Icdx  is  denoted  "by  a. 


VctM  p 

'  -A* 

2)B3c 


(12) 


Note . that  even  though  the  expression  for  0Na  in  equation  (12)  does 

not  contain  the  integrated  effects  of  higher-order  terms  in  co,  it  remain 
sinusoidal  vith  respect  to  time. 
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In  order  to  obtain  it  is  only  necessary  to  change  the  definition 

of  W(|,q)  in  equation  (10).  In  this  case 


W(  |,ii )  =  io»0(  i  -  Xq ) 


(13) 


Substituting  the  expressions  from  equations  (9)  and  (13)  into  equation  (8) 
and  letting  z  approach  zero  yields 


I  _  i<m 


x  pn2  ( 1  “  xo)e  cP 


ioM, 

— (X"?) 


Pc 


COS 


vt”  -  "1 jfs  -  ”) 


i’ll 


\f( "  -  \)("2  - ") 


dq  d|  (11+) 


which,  after  expansion  to  first  power  of  cd,  yields 


0Na  ~ 


ion- 

*P 


"12  (l-Vl1-iT 

c(3 


E 


(x-  |) 


dr)  d| 


or 


0Na  =  ~  x(x  -  2xQ) 


=  ^  (x2  -  2xxq) 


(15) 


Equations  (12)  and  (15)  may  "be  substituted  into  equation  (5)  to  obtain  the 
complete  expression  for  0^  at  the  upper  surface  of  the  wing. 

Derivation  of  0^ • —  As  pointed  out  in  the  introduction^  distributions 

of  doublets  may  be  used  to  satisfy  the  boundary— value  problem  in  the  region  of 
mixed  supersonic  flow  (region  T).  The  type  of  doublet  required  is  that  with 
axis  normal  to  the  reference  plane.  It  will  be  recalled  that  the  potential 
of  such  a  doublet  may  be  obtained  from  the  potential  of  a  source ,  located 
in  the  reference  plane ,  by  a  partial  differentiation  of  the  source 
potential  with  respect  to  the  direction  normal  to  the  reference  plane. 
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If  such  a  differentiation  is  applied  to  the  potential  of  a  moving  source 
in  the  plane  (equation  (7)),  the  following  expression  that  may  he  shown 
to  satisfy  equation  (l)  is  obtained: 


where  now  W(£,q)  refers  to  doublet  strength. 


In  order  to  obtain  0rpa  for  the  oscillating  flat  wing  it  is  convenient 
to  assign  to  the  functions  w(t  -  r ±)  +  w(t  -  tp)  and  W(|,t))  the 
expressions  given  in  equation  (9)  and  equation  (10),  respectively. 
Substituting  these  expressions  into  equation  (l6)  gives 


ft 


ir 
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Expanding  the  terms  involving  cn  in  equation  (17)  into  a  power  series 
to  the  first  order  of  co  gives 


02  ft:  — 2VaP^z 


n-i) 


("  -  \)(\  -  "j 


=  -2VP2 


l"  ~  "iX 


•n2  —  ”n 


3/2  cp: 


j(n  -iiX’ 


Each  of  the  two  terms  within  the  traces  in  equation  (l8)  can  he  shown  to 
satisfy  equation  (l)  to  the  first  order  of  go.  Furthermore  it  can  he  seen 
that,  to  the  first  order  of  the  potential  of  a  moving  doublet  may  he 
considered  as  a  superposition  of  two  doublets  of  the  stationary  type;  one 


of  strength  proportional  to  a 


and  the  other  of  strength  proportional 


to 

P 

cp 


The  potential  0^a  composed  of  distributions  of  doublets  of  the  types 
in  equation  (l8)  may  be  written 


(19) 
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where  the  integration  is  to  he  extended  over  the  wing  surface  S  included 
in  the  fore  cons  with  vertex  at  (x,y,z)  (see  fig.  2(b)  for  the  limiting 
case  z  =  0)  and  where  and  F2  are  functions  describing  the  manner  in 

which  the  doublets  are  to  be  distributed  in  the  region  S  and  must  be  deter¬ 
mined  in  such  a  manner  that  the  boundary  conditions  (a)  and  (b)  will  be 
satisfied. 

The  normal  procedure  of  determining  Fj_  and  F2  would  be  to  formulate 

an  integral  equation  by  imposing  the  boundary  conditions  on  equation  (19) ; 
however,  this  procedure  is  tedious  and  difficult.  It  can  be  circumvented 
here  by  making  use  of  known  solutions  to  similar  problems.  The  integrand 
of  the  first  integral  in  equation  (19)  is  recognized  to  be  the  same  as  that 
in  the  integral  form  of  the  velocity  potential  for  a  similar  wing  at  constant 
angle  of  attack,  and  the  integrand  of  the  second  integral  is  recognized  to 
be  the  same  as  that  in  the  integral  form  of  the  velocity  potential  for  a 
similar  wing  undergoing  a  constant  rotation  about  its  leading  edge. 
Furthermore,  the  boundary  conditions  for  these  cases  can  be  stated  precisely 
as  conditions  (a)  and  (b),  stated  herein,  for  the  oscillating  wing.  The 
distribution  function  and  the  velocity  potential  at  the  wing  surface  for 
the  rectangular  wing  at  constant  angle  of  attack  can  be  readily  obtained 
from  the  known  distribution  of  lift  for  this  case.  The  distribution 
function  and  velocity  potential  at  the  wing  surface  for  the  rectangular 
wing  undergoing  a  constant  rotation  about  its  leading  edge  can  be  obtained 
from  the  wing  at  constant  angle  of  attack  by  linear  superposition.  The 
distribution  functions  and  velocity  potentials  at  the  wing  surface  for 

both  these  cases  are  derived  in  the  appendix.  The  value  of  0 

”Ta 

(equation  19)  at  the  wing  surface  can  be  easily  deduced  from  these  results; 
that  is. 


&  fa 

P* 


l/Py(x  -  Py)  +  x  sin-1 


VMa 

cp3it 


vV(x  -  Py)  +  x2 
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2Va 


\j  Py(x  -  3y)  +  x  sin- 


r-1 


M2a 

3«33 


(x  +  2py)  \/  Py(*  -  Py) 


+  3x2sin  1 


(20) 


In  order  to  obtain  0T(^  it  is  convenient  to  define  V(  ^t)  )  as  in 
equation  (13).  Substituting  equations  (9)  and  (13)  into  equation  (16) 
yields 


02  =  “ 

oz 


■2a. (|  -  x0) 


-i^(x-|) 

*  cp^  cos 


a> 

Pc 


'/(l  -  ^(ig  -  ’i) 


\j("  ^  ") 


(21) 


Carrying  out  the  indicated  differentiation  and  expanding  to  the  first 
power  of  m  results  in 


02  =  -2p2za  ^  ^ 


[("  -  \)(\  -  ') 


3/2 


2|32zd< 


•0 


V. 


3/2 


3/2 


(22) 


The  two  expressions  in  the  right  side  of  equation  (22)  have  the  same  form 
as  the  expressions  in  the  right  side  of  equation  (l8).  The  following 
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integral  expression  for  0^  is  therefore  of  the  same  form  as  the 
integral  expression  for  0Ta  (equation  19)  and  the  functions  F^ 

and  can  he  determined  by  the  procedure  given  for  determining 

and  F2: 


The  expression  for  0rp^  at  the  wing  surface  can  therefore  be  deduced  from 
the  results  in  the  appendix  to  be 


(23) 

Equations  (20)  and  (23)  may  be  substituted  into  equation  (6)  to  obtain  0-p 
at  the  wing  surface. 

Aerodynamic  moment.—  The  local  pressure  difference  between  the  upper 
and  lower  surfaces  may  be  written  as 
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IT 


The  aerodynamic  moment  on  the  entire  wing  about  the  axis  x  =  Xq,  measured 
in  the  positive  clockwise  direction,  is  then  (see  fig.  1) 


Ma  =  2 


(x  -  dy  dx 


+  2 


//< 


-  X(^  Pp  dy  dx  (25) 


If  all  linear  dimensions  involved  in  equations  (2k)  and  (25)  are  referred 
to  the  chord  2b,  these  equations  become,  in  a  sense,  nondimensional  and 
can  be  written 


(27) 


Putting  in  the  appropriate  limits  of  integration  (see  fig.  2)  results  in 
the  following  equation  for  the  total  aerodynamic  moment : 


M  =  16b3 

OC/ 


>h 


10  Jx/P 


(x  -  xQ)  Pjj  .dy  dx  +  16b3 


(* 


(28) 


If  the  nondimensional  forms  of  p^.  and  are  computed  from 

equations  (12),  (15),  (20),  and  (23)  and  substituted  into  equation  (28), 
there  is  obtained  a  complex  expression  for  that  may  be  written  as 


=  — Spb-^^k  ocA^M^  +  iM^ 


(29) 
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where 


U  =  —1— 
3  302k2 


k  "  2^)  "  A  (2  -  3*°) 


(30) 


is  the  in-phase  component  of  the  aerodynamic-moment  coefficient,  and 


%  =  -fr 

3^k 


2P3(i  -  3xo  +  3xq2)  -  0(2  -  3*o)  -  2p2xo^3xn — !l Jtfco. 3) 

2A 


(31) 


is  the  out— of —phase  or  torsional-damping  aerodynami c-moment  coefficient. 

The  expression  for  in  equation  (31)  represents  the  central  result 

of  the  present  investigation  “because  the  conditions  of  torsional  stability 
or  instability  depend  directly  on  the  sign  and  magnitude  of  M^.  It  will 

be  noted ,  in  equation  (31)*  that  may  undergo  changes  in  sign  as  any 

one  of  the  parameters  P,  xQ,  and  A  is  varied. 

Although  equation  (31)  has  been  derived  for  the  condition  that  the 
intersection  of  the  Mach  cones  from  the  foremost  points  of  the  side  edges 
is  off  the  wing,  it  can  be  shown  that  it  remains  valid  as  long  as  the  Mach 
cone  from  one  side  edge  does  not  intersect  the  opposite  side  edge.  When  a 
side  edge  is  intersected  by  the  Mach  cone  from  the  opposite  edge  different 
velocity  potentials  from  those  explicitly  derived  herein  are  required  for 
portions  of  the  wing  behind  the  points  where  the  Mach  cone  intersects  the 
edge.  The  derivation  of  these  potentials  is  not  considered  here. 

Total  torsional— damping  coefficient  and  some  selected  calculations  of 
stability  conditions.—  If  mechanical  damping  is  assumed  to  exist  about  the 
axis  of  rotation,  for  example  as  the  hinge  friction  in  an  aileron  instal¬ 
lation,  and  if  this  damping  is  converted  to  coefficient  form  as  in 
references  3  and  9  and  combined  with  the  aerodynamic  damping,  the  total 
torsional-damping  coefficient  may  be  written  as 

V  =  M4(0,Xo,k,A)  +  nra2(^2ga 
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or 


Mi/ 


(32) 


Negative  values  of  indicate  dynamic  instability  that  may 

correspond  to  single— degree  torsional  flutter.  Positive  values  of  M/ 

indicate  stable  conditions  and  between  the  stable  and  unstable  conditions, 
that  is  when  M^?  vanishes,  a  borderline  state  of  unstable  equilibrium 

separating  damped  and  undamped  torsional  oscillations  occurs.  Since  the 
term  involving  in  equation  (32)  is  always  positive,  negative  values 

of  M^?  and  the  vanishing  of  M^1  are  associated  with  negative  values 

only  of  (equation  ( 31) )  - 

With  the  exception  of  aspect  ratio  A,  the  effects  of  the  individual 
parameters  in  equation  (32)  are  presented  in  references  3  and  4.  The 
effect  of  aspect  ratio  for  some  particular  values  of  the  other  parameters 
is  briefly  discussed  in  the  following  paragraphs. 

The  range  of  values  of  Xq  and  0^  (and  M^)  for  which  M/  vanishes 

when  ga  is  assumed  to  be  zero  and  A  has  different  selected  values  is 

shown  plotted  in  figure  3-  The  regions  inside  the  different  curves  give 
the  range  of  values  of  Xq,  0^,  and  in  which  instability  might  occur. 

The  curve  for  infinite  aspect  ratio  agrees  with  that  given  in  figure  21  of 
reference  3  and  that  given  in  figure  2  of  reference  4.  It  will  be  noted 
that  decreasing  the  aspect  ratio  has  a  highly  stabilizing  effect. 

The  dashed  portions  of  the  curves  in  figure  3  correspond  to  aspect 
ratio  and  Mach  number  combinations  that  cause  the  Mach  cone  from  the 
foremost  point  of  one  tip  to  intersect  the  opposite  tip  ahead  of  the 
trailing  edge. 

In  figure  4  the  quantity  M^  multiplied,  for  convenience,  by  33^k  is 

plotted  as  a  function  of  M^  for  the  values  of  A  selected  in  figure  3 • 

The  ordinates  of  these  curves  are  proportional  to  the  amount  of  negative 
aerodynamic  damping  available  for  given  Mach  numbers  and  for  given  values 
of  the  frequency  parameter  k.  .  The  stabilizing  effect  of  decreasing  aspect 
ratio  is  apparent  in  this  figure. 

Corresponding  values  of  the  flutter-speed  coefficient  V/boc^  and 
aspect  ratio  A  for  which  M^1  vanishes  for  some  selected  values  of  the 

parameters  xQ,  g^,  M2,  p,  p  and  ra2  are  plotted  in  figure  5. 
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The  regions  above  these  curves  give  the  range  of  values  in  which  insta¬ 
bilities  might  occur. 

It  should  be  recalled  that  thickness  effects  have  not  been  considered 
in  the  present  discussion.  Second— order  thickness  effects  may  be  of 
considerable  significance  in  regard  to  undamped  torsional  oscillations  of 
wings  of  finite  or  infinite  aspect  ratio. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 

Langley  Air  Force  Base,  Va.,  March  28,  1949 
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APPENDIX 

VELOCITY  POTENTIALS  FOR  TIP  REGION  OF  RECTANGULAR 
WING  IN  STEADY  SUPERSONIC  FLOW 


It  can  be  shown,  as  in  reference  10,  that  at  points  in  the  xy- plane 
the  perturbation— 'velocity  potential  in  a  supersonic  stream  due  to  doublets 
in  the  xy— plane  distributed  according  to  a  distribution  function  f(|,Tj) 
is  proportional  to  the  distribution  function;  that  is. 


w/  Vapz  ff  f(M)dT|  d| 

0(xjYjO)  -  H®1  n  II  “lo/o 

Jw  |u-  D2-U2(y-n)2-3^j3/ 

=  ±.itf  (x,y)  (Al) 


where  the  positive  sign  applies  to  the  upper  surface  of  the  region  S  and 
the  negative  sign  applies  to  the  lower  surface;  that  is*  the  direction  from 
which  z  0.  When  the  appropriate  distribution  function  can  be  found,  such 
doublet  distributions  can  be  used  to  satisfy  boundary-value  problems,  in  the 
xy-plane,  for  vanishingly  thin  airfoils  at  vanishingly  small  angles  of 
attack.  In  the  case  of  a  rectangular  wing  placed  in  a  supersonic  stream 
flowing  from  the  negative  x-direction,  with  its  leading  edge  normal  to  the 
free— stream  direction,  the  added  velocity  in  the  free— stream  direction  in 
the  tip  region  can  be  obtained  from  the  expression  derived  for  the  added 
pressure  difference,  in  this  region,  in  reference  11.  In  the  present 
notation  this  expression  for  pressure  difference  is 


2pT  M ,  sm-1  ./n 

Sx  Pit  »  x 


or 


(A2) 
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If  0  ia  eliminated  from  equations  (Al)  and  (A2),  the  following  differ¬ 
ential  equation  for  f (x,y)  is  obtained: 


jt 


8f 


(A3) 


The  expression  for  rtf(x,y)  or  JZS(x,y)  can  now  he  obtained  by  a  partial 
integration,  with  respect  to  x,  of  equation  (A3).  This  integration  gives 


KfU^y)  =  0(x,y) 


—  Py) +  x  sin 


-1 


(Ah) 


For  the  rectangular  wing  rotating  about  its  leading  edge  the 
distribution  function  for  the  tip  region  can  be  obtained  by  a  linear 
superposition  of  the  function  given  in  (Ah);  that  is. 


—it 


f(x  -  |)d|-  = 


i)  -  PyJ  +  (x  -  1)  sin 


Va 

P* 


-  2Py 
3 


\/py(x  -  j3y)  +  x2  sin  1 


(A5) 


If  m  equation  (Ah)  a  is  replaced  by  a - J  and  if  in  equation  (A5) 

.16  cP 

ct,  is  replaced  by  — the  sum  of  the  resulting  equations  will  yield 

cp 

equation  (20)  of '  the  text.  Similarly,  if  in  equation  (Ah)  Va  is 
replaced  by  — XqO,  and  if  in  equation  (A5)  Va  is  replaced  by  a,  the 

sum  of  the  resulting  equations  will  yield  equation  (23)  of  the  text. 
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(b)  Section  y  =  y^  (xz-plane). 

Figure  1.-  Sketch  Illustrating  chosen  coordinate  system 
of  airfoil  oscillating  in  torsion. 
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(b)  Mixed  supersonic  region. 

Figure  2.-  Sketch  illustrating  areas  of  integration  for  "purely  supersonic 
and  "mixed  supersonic"  regions  of  flow. 


Stable  region 
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damping  vanishes  for  some  selected  values  of 
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Figure  4.-  Curves  showing  the  variation  of 
negative  aerodynamic  torsional  damping 
with  M2  for  selected  values  of  A 
and  Xq  • 


j 


(3)3  k)M4 
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Figure  4.-  Continued. 
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